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Abstract 

We discuss Mott insulating and metallic phases of a model with e g orbital degeneracy to 
understand physics of Mn perovskite compounds. Quantum Monte Carlo and Lanczos 
diagonalization results are discussed in this model. To reproduce experimental results on 
charge gap and Jahn- Teller distortions, we show that a synergy between the strong corre- 
lation effects and the Jahn- Teller coupling is important. The incoherent charge dynamics 
and strong charge fluctuations are characteristic of the metallic phase accompanied with 
critical enhancement of short-ranged orbital correlation near the insulator. 

1 Introduction 

Transition metal oxides with perovskite structure offer many prototypical systems which 
show various types of strong correlation effects [[j]]. Mott insulating states are widely 
observed phenomena in these compounds at integer fillings of 3d electrons. A typical 
property of metals found near this Mott insulating state is that they are "bad metals" 
where transport and optical measurements show strong dampings and dominance of in- 
coherence in the charge dynamics. The optical conductivity generically has a broad and 
long tail up to the order of eV, which is in contrast with the prediction of the simple 
Drude theory. The available photoemission data also in general indicate strong damping 
effects and smallness of quasi-particles weight near the Fermi surface. A common origin 
of this incoherence may be proximity of the Mott insulator with dominance of incoherent 
excitations generated by disorder of spin or orbital degrees of freedom at low tempera- 
tures. In other words, spin and orbital fluctuations scatter charge excitations incoherently 
if they are not ordered. Among the transition metal oxides, Mn oxides, Lai-^Sr^MnOs 
have a unique character, where ferromagnetic metals exist with a perfect spin polariza- 
tion at low temperatures due to the double exchange mechanism. Even under the spin 
polarization, the Mn compounds show a similar or even more incoherent charge dynamics 
than other transition metal oxides 0, |3], [|, [J . This implies importance of orbital fluctu- 
ations presumably coupled with lattice fluctuations, because the spin degrees of freedom 
are not able to contribute to the incoherence. The orbital degrees of freedom coupled to 
lattice are also expected to play an important role in the insulating phase through the 
cooperative Jahn- Teller ordering. 

An important and puzzling feature in the ferromagnetic metallic phase of Lai-^Sr^MnOa 
is that the Drude weight relative to the noninteracting system is very small (~ 0.01-0.05) 
while the specific heat coefficient 7 is also small (~ 2-5 mJ/moleK 2 ) ||. The former is a 
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typical indication of the incoherent charge dynamics, which is likely to be a consequence of 
orbital and lattice fluctuations. Naively, these fluctuations and resultant residual entropy 
are expected to increase 7, which is not the case in this compound. 

In this paper, using numerical methods, we analyze roles of orbital degrees of freedom 
in a spin polarized two-dimensional model J7L Rl . To understand fluctuation effects in de- 
tail we employ a quantum Monte Carlo (QMC) and Lanczos diagonalization methods for 
numerical calculations. Our numerical results are summarized as follows: At half filling 
where the e g orbital is occupied by one electron per site on average, interorbital on-site 
Coulomb repulsions U opens the Mott charge gap . The opening of the gap is, how- 
ever, substantially slower than the Hartree-Fock prediction of the present model and also 
than the QMC results of the corresponding single-band Hubbard model. The staggered 
orbital order appears at the Mott insulating phase, while the order is destroyed quickly 
if the system is metallized by doping of carriers. The short-ranged orbital correlation 
shows critical enhancement when the Mott insulator is approached. This is similar to 
the criticality observed in the single-band Hubbard model, where spin correlations are 
enhanced instead of the orbital correlations. In accordance with critical fluctuations of 
the orbital correlations, the Drude weight is more strongly suppressed with decreasing 
doping concentration than the total intraband weight. 

When the coupling of the orbital to Jahn- Teller distortions is introduced, the charge 
gap is strongly enhanced through the optimized Jahn- Teller distortion. We optimize 
the static Jahn- Teller distortion to lower the total energy including electronic, lattice 
and their coupled part. This corresponds to a classical treatment of the lattice. A 
reallistic choice of electron-lattice coupling suggests that the Jahn- Teller distortion and 
the on-site Coulomb repulsion are both important for the formation of the charge gap 
in LaMnOs through synergetic effects of these two. With this choice of parameters, the 
charge gap amplitude, stabilization energy and amplitude of the Jahn- Teller distortion are 
all rather quantitatively reproduced. In the metallic phase, we argue that the Jahn- Teller 
coupling induces an instability to phase separation near the Mott insulator. We also 
discuss effects of three dimensionality, which retains small Drude weight but contributes 
to the suppression of the specific heat coefficient 7 through an anisotropic scaling [[J. 

2 Model 

We study an orbitally degenerate model of the manganese oxides with electronic, lattice 
and their coupled terms described by 

H — H e l + Hph + Hel-ph- (1) 

The first term 7i e \ represents the electronic part which contains the kinetic energy and 
the Coulomb interaction. A Mn 3+ ion in LaMnOs has four 3d electrons. Three electrons 
occupy t 2g orbital while the last electron occupies one of the doubly degenerate e g orbital, 
where all of the four electrons have parallel alignment of the spins due to strong Hund's 
rule coupling. In the ferromagnetic and metallic phase, the spins are perfectly polarized 
with the moment ~ 3.8/^b at low temperatures and the spin degrees of freedom do not 
play a role. Therefore a model which ignores the spin degrees of freedom would be a 
good starting point to discuss low-energy physics of the Mn compounds in the metallic 
phase. Moreover, the Mott insulating phase has A-type antiferromagnetic order in which 
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the spins are ordered in parallel in a two-dimensional (2D) plane. Therefore, a 2D plane 
isolated from the perovskite structure always has a ferromagnetic polarization near the 
Mott insulator. We will discuss later that the 3D coupling only contributes to a detailed 
low-energy structure in the insulating phase, as naturally expected from magnetic as well 
as orbital orderings with strong 2D anisotropy. Although the orbital long-ranged order is 
lost, the strong 2D anisotropy may also be retained in metals, and we study the metal- 
insulator transition in the 2D model as a starting point. We later discuss 3D effects on the 
metal-insulator transition based on these 2D results. We also note that our model may 
have some relevance to layered Mn compounds [H|, which consist of stacking of single or 
double Mn0 2 layers with rather small crystal- field splitting of d x 2_ y 2 and d^ z 2_ T 2. All of 
the above problems may be studied in an isolated 2D plane at least as a good starting 
point. We employ such an appropriate 2D model by ignoring the spin degrees of freedom. 
In this model, we still have orbital degrees of freedom for e g electrons. Therefore the 
electronic part of the Hamiltonian Ti, e \ describes e g electrons with double degeneracy as 



H e i = - U,j,uy{4 u c ju > + H.c.) 

+U^2{na - tjX^- ( 2 ) 

i iv 

where the creation (annihilation) of the single-band electron at site % with orbital v is 
denoted by c iv {c iv ) with n iv being the number operator n iu = c\ v Ci U . The Hamiltonian 
looks similar to the ordinary Hubbard model if the orbitals are replaced with spins. 
However, an important difference from the ordinary Hubbard model is that the transfer 
has an anisotropy and off-diagonal elements with dependence on the d x 2_ y 2 and d^ z 2_ r 2 
orbitals as well as on spatial directions, where we label the orbital d x 2_ y 2 and d% z 2_ T 2 as 
v — 1 and 2, respectively below. In 2D configuration, the nearest-neighbor transfer is 
scaled by a single parameter t as t u = |i , t 22 = \to } ti2 = hi = ±-^0, where, in ±, 
+ (— ) is for the transfer to the y(x) direction. We note that the interaction U is not a bare 
interorbital Coulomb repulsion £7 12 but £/ 12 — J12, where J i2 is the Hund's rule coupling 
between orbitals v = 1 and 2. In this paper we take the energy unit as to- 

In the absence of U, the diagonalized two bands have dispersions both with the band- 
width 4to, where the energy of one band is 2to higher as a whole than the other. We have 
to note that the perfect nesting condition is satisfied at half filling in the noninteracting 
isolated 2D plane although it is not in the real 3D structure. In the 2D lattice at half 
filling, the lower band is at three-quarters filling while the upper band is at a quarter 
filling, both of which have square shapes of the Fermi surfaces. The nesting vector (it, it) 
connects two Fermi surfaces belonging to different bands. Because of this nesting con- 
dition in the 2D lattice, it may have stronger tendency to the Mott insulator with the 
orbital order at half filling than the 3D structure. 

The electronic degrees of freedom are coupled with lattice through several different 



types of distortions | IT]| . We consider the Jahn- Teller coupling in the Hamiltonian (|I|), be- 
cause it linearly and strongly couples with the e g orbital fluctuations and ordering, which 
is the main subject of this paper. The Jahn- Teller distortion coupled to two e g orbitals is 
represented by two orthogonal modes of oxygen displacements «j at the octahedron corner 
surrounding a Mn ion at the i-th site. A linear combination of these two modes Q 2 and 
Q3 may represent any Jahn- Teller-type distortion coupled to the e g orbitals, expanded as 
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a two-component vector as Vj = — 2wj (cos 2^, sin2#j), where 9i determines the angle of 
the linear combination. The Q2 and Q3 modes linearly couple with the orbital polariza- 
tion Tf and Tf, respectively, where the pseudospin operator T/*(/i = x, y, z) is defined by 
Tf 1 = \ &w' c iv c iu', with a being the Pauli matrix. Then the electron-lattice coupling 
and the lattice elastic energy are given as 

Ki-ph = gi ■ I;, 7-iph = kJ2 u li (3) 

i i 

where g« = gvj and Ij = (T?,T?). The electron-lattice coupling constant and the spring 
constant for the Jahn- Teller lattice distortion are represented by g and k, respectively. 
We treat Ui and 9i as a classical variable [|l2j. In the following the possibility of static and 



uniform Jahn- Teller distortions u = Ui with experimentally-observed staggered pattern 
Oi = ~(— l)l r "l + l ri al is considered by minimizing the total energy. 

When the Jahn- Teller distortion is nonzero, the transfer amplitude to is also modified 
through the oxygen distortion. Empirically it is expected that the Mn-0 transfer changes 
by a factor (1 ± 2-Uj)~ 7 / 2 for the oxygen distortion U{ |13[]. In this paper, we take the 



length unit as the Mn-Mn distance. Then, although this has minor effects on our results 
shown below, we take the modified Mn-Mn transfer as to = to(l — 4'U 2 ) -7 / 2 . 

The model ([I]) is suited to understand two different problems. The first is the basic 
electronic structure where orbital-lattice coupled systems primarily determines the insu- 
lating charge gap while the intersite spin correlations play secondary roles. The second 
is the incoherent charge dynamics in metals at low temperatures where spin degrees of 
freedom are indeed absent due to the perfect polarization. 

We employ the projector quantum Monte Carlo (PQMC) [|i~4| , |15f and the Lanczos 
diagonalization (LD) methods for the model ([]])- (H) to calculate physical properties in 
the ground state. In the PQMC, the minus sign problem poses a restriction to tractable 
range of the parameters. We show the cases in which convergence to the ground state is 
reached without this problem and the ground state properties are obtained with needed 
accuracy. In the next section, the charge gap and orbital order and Jahn- Teller distortion 
are presented after taking the extrapolation to the infinite system size (thermodynamic 
limit). In the LD, system size is limited by the computer memory size. We calculate two 
sizes, 4 by 4 and v^lO by y/lO lattices. The boundary condition is chosen so as to realize 
the ground state, where phase difference is varied in x and y directions. Note that = 
and = 7r correspond to the periodic and antiperiodic boundary conditions, respectively. 

As realistic choices of parameters in the Hamiltonian, we take the following 0, [16], [H], 
|l8i : The effective Mn-Mn transfer to is estimated around 0.5eV while the on-site Coulomb 
repulsion is varied in the range of < U/to < 16 to see the U dependence. The elastic 
constant k is roughly estimated from the Mn-0 stretching phonon frequency (~ 70 meV) 
as several tens eV in the length unit of Mn-Mn distance 0, |3|, [L2| . In our calculation, we 
take k/to = 100. In the same length unit, the Jahn- Teller coupling g is taken as g/t = 10, 
which can be roughly estimated from the electrostatic energy between an 2 ~ ion and an 
electron on different types of 3c? orbital on the neighboring Mn site. 

3 Results 

The results are discussed in the following order: First we discuss results at half filling. 
The ground state properties of Ti e \ without the Jahn- Teller coupling are presented. We 
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next show how the results are modified when we introduce the Jahn- Teller coupling. In 
the second part, carriers are doped and the metallic phase is studied. Away from half 
filling, we mainly show properties in the absence of the Jahn- Teller coupling, although its 
effects are briefly discussed in the next section. 

At half filling without the Jahn- Teller coupling, the charge gap and the orbital order 
were calculated by the PQMC method. As we see in Fig. [I], the charge gap amplitude A c 
is substantially smaller than the Hartree-Fock prediction. Although the nesting condition 
is satisfied in the 2D model, the growth of the charge gap with increasing U is slow. This 
slow growth is also confirmed from the comparison with the PQMC results of the ordinary 
single-band Hubbard model. The charge gap at U/to = 4 is around 0.66 for the single- 
band Hubbard model while it is around 0.1, in the present case. To reproduce the Mott 



charge gap in the experimental result of LaMnOa, A c ~ leV [111, we need much larger U. 
This may also be true in real 3D structure because the Mott gap opening would be even 
slower due to the absence of the perfect nesting. The slow growth of the gap may be due 
to fluctuations and competition between interband off-diagonal and intraband hoppings 
which both satisfy the perfect nesting condition while the orbital order is stabilized only 
through the interband part. The orbital order is also plotted in Fig. [I], where the order 
grows quickly beyond some threshold U/to ~ 4. We note that this orbital polarization 
has a uniaxial symmetry. 

When the Jahn- Teller coupling is introduced with static and uniform Jahn- Teller dis- 
tortion u with the staggered pattern 6i being allowed, the ground state is achieved at a 
nonzero u. Figure || shows the PQMC results of the ground state energy at U/to — 4 
as function of static and uniform distortion u. With this choice of parameters k and g, 
the U dependence of the optimized distortion is plotted in Fig. |3|. The experimentally 
observed value u = 0.035 for LaMnOs is close to the values of the optimized distortions 
plotted in this figure. Note that the optimized distortion u at U = is nearly zero. The 
small deviation between our numerical results around U — 4 and the experimental results 
would possibly be due to the uncertainty in the choice of k and g. We will discuss this 
point later. The U dependence of the charge gap at the optimized distortion is plotted 
in Fig. |3|. As compared to the case g = 0, opening of the charge gap is much quicker 
with increasing U. We note that the gap A c at U = is around 0.02t and far below 
the experimentally observed value ~ 2t . The quick opening of the gap with increasing 
U comes from synergetic and nonlinear effects between electron correlations U and the 
Jahn- Teller distortions. The Jahn- Teller distortions and the on-site Coulomb repulsion 
cooperatively stabilize and enhance the charge gap of the insulating state through the dy- 
namic coupling of orbital and lattice polarizations. The experimentally observed charge 
gap A c ~ leV (A c /i ~ 2) is numerically reproduced if we take U/to at around 4 and 5. 

The ground state energy difference, £jt, between the case with optimized Jahn- Teller 
distortion and that at u = roughly measures the stabilization energy of the Jahn- Teller 
distortion. This is around 0.08eV (~ 0.1-0.2to) if we take U/to — 4-5 as seen in Fig. It 
is again favorably compared with the Jahn- Teller transition temperature ~ 800K. 

From these analyses, a consistent picture to explain the experimental and the present 
numerical results emerges : (1) The charge gap of the order of leV is determined by a syn- 
ergetic effect of the on-site Coulomb repulsion and the atomic Jahn- Teller distortion. (2) 
The Jahn- Teller transition temperature accompanied by the orbital long-ranged order is 
roughly determined from intersite coupling for combined Jahn- Teller and orbital exchange. 
The intersite orbital correlation and intersite correlations of the Jahn- Teller distortion are 
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both developed in the energy scale of ~ O.leV. (3) The interlayer antiferromagnetic spin 
correlation grows at much lower energy scale around 0.02 eV. 

In our analyses, although we have fully taken into account the orbital correlations 
including short-ranged and dynamical fluctuations, we have not considered dynamical 
effects for the Jahn- Teller distortions. This will have the following effects: When a hole 
is doped, the Jahn- Teller distortion may relax around the hole site to form a polaron. 
The one-hole energy then decreases and the charge gap A c may thus to some extent be 
reduced from the present value. Next, because Ej T is the energy difference between the 
Jahn- Teller distorted and undistorted phases, it may overestimate the energy scale of the 
Jahn- Teller transition temperature Tjt because £"jt contains not only the intersite Jahn- 
Teller ordering energy but also the on-site Jahn- Teller polarization energy, while Tjt is 
determined basically from the intersite part [TTJ]. Presumably, larger values of k and g 
than the present choices would be required to reproduce A c and T JT if we consider the 
short-ranged Jahn- Teller correlation effects. We need further studies on this short-ranged- 
correlation effects and also on the uncertainties in the choice of the parameter values for 
k and g. 

We next discuss the doped cases. When the Jahn- Teller coupling is absent, the orbital 
order clearly observed in the half-filled Mott insulator is quickly destroyed with increasing 
doping concentration 5. In Fig. ^, equal-time structure factor of the orbital correlation 
T X (Q) is plotted at Q = (tt, tt) as function of hole concentration 5. Here we define 
T M (k) = iVg" 1 YsijiTt^Tj 1 ) exp (ikr^) (/i = x,y,z) for an iV s -site system. This indicates 
that T X (Q) is scaled by 5" 1 hence the orbital correlation length is scaled by the mean 
hole distance <5 -1 / 2 . This is the same scaling as the antiferromagnetic correlation in the 2D 
single-band Hubbard model [20]. The doping concentration dependence of the chemical 
potential is plotted in Fig. [j| This implies that the chemical potential is proportional to 
5 2 , which results in the scaling of the charge susceptibility Xc = dn/d/j, oc This is 
again the same scaling as the 2D single-band Hubbard model |20| , |2l| . 

These suggest that the critical enhancement of the orbital correlation plays a simi- 
lar role to spin correlations in the single-band Hubbard model and the metal-insulator 
transition is controlled by a similar criticality if the Jahn- Teller coupling is switched off. 
Although it is not numerically confirmed, it is tempting to ascribe the enhancement of 
the charge susceptibility to flat dispersions generated at some particular points of the 
momentum space in the single-particle excitation near the Fermi level at small 5, as is 
indeed observed in the single-band Hubbard model [p2 |. This may generate dynamical 



fluctuations for the phase separation or dynamical charge-order fluctuations related to 
the location of the flat dispersions. 

The optical conductivity calculated by the LD method shows reduction of the Drude 
weight with decreasing 5. This reduced Drude response coexists with a conspicuous 
intraband incoherent response at finite frequencies and also transition to the upper (or 
lower) Hubbard band at high frequencies as seen in Fig. |[ In this figure, U is taken 
larger than the realistic values to separate the intraband incoherent response from the 
inter-Hubbard-band transition. In the experimental situation, these two finite-frequency 
weights are merged to a single structure. A large incoherent response at finite u within the 
Mott gap may be a consequence of carrier scattering by short-ranged orbital correlations. 
The optical conductivity will be studied elsewhere in more detail [§]. 
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4 Discussion 



Quantum Monte Carlo results indicate that for realistic choice of the parameters of 
LaMn03, the insulating charge gap (~ leV) at half filling is a consequence of synergetic 
and cooperative enhancement by the "Mott gap" and "Jahn- Teller gap". The energy 
scale of intersite correlations for the Jahn- Teller distortions and the orbital correlations 
is one order of magnitude smaller than the charge gap energy scale (~ leV) . The energy 
scales of antiferromagnetic correlations and the three-dimensional effects are two orders of 
magnitude smaller than the charge gap scale. The experimental results for the charge gap 
A c and the Jahn- Teller distortion u are favorably compared with our results at U/to = 4- 
5,g/to ~ 10, k/t ~ 100 and t ~ 0.5eV. To make more precise comparisons possible 
between numerical and experimental values, a more scrutinized tuning of the parameters 
would be required. In the absence of the Jahn- Teller coupling, the orbital long-ranged 
order stabilized at half filling is destroyed quickly upon doping while the short-ranged 
orbital correlation is critically enhanced with 5 — > in metals. The charge susceptibility 
shows critical enhancement with 5 — > implying the formation of a flat dispersion at 
some particular points in the momentum space of the single-particle excitation. 

When the Jahn- Teller coupling is introduced for the doped cases, our numerical results 
suggest that, in the absence of long-ranged Coulomb interaction, the phase separation may 
take place in a finite range of doping concentration typically between 5 = and 5 ~ 0.2. 
This phase separation is resulted from a combined effect of the already enhanced charge 
susceptibility at g = and the Jahn- Teller effect. Note that the phase separation takes 
place if the charge susceptibility diverges. A similar result of the phase separation was 



suggested in a study of a one-dimensional model p3| . With longer-ranged Coulomb inter- 



action, charge ordering may take place instead of real phase separation. The Jahn- Teller 
coupling promotes such tendency toward real-space structure in non-half-filled systems. 
This problem will be discussed elsewhere in detail [24j| . 

If a dynamical charge fluctuation and flat dispersions appear, strong reduction of the 
Drude weight is a natural consequence as already clarified in the studies of the 2D t-J 
model p5" |. The scaling theory of 2D system also predicts critical enhancement of the 



specific heat coefficient 7 []26|] . This is however, experimentally not observed as discussed 
in §1. This problem was examined in terms of three-dimensional effects When a flat 
dispersion in the form of ~ k A appears in the intraplane direction while the dispersion 
is normal in the interlayer direction, the specific heat coefficient follows 7 oc 5° and does 
not show critical enhancement while the Drude weight follows D oc <5 4 / 3 . This seems 
to be consistent with the experimental observation. We need further studies to fully 
understand the microscopic basis in detail. The momentum dependence of the single- 
particle excitations is an important issue to be explored. 
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Figure 1: QMC results (circles) without the Jahn- Teller coupling for U dependences of 
(a) the charge gap and (b) the staggered orbital polarization. The data are shown in the 
limit of iVg — > oo, which are obtained from size extrapolations of the data on the sizes 
iVg = 4 x 4 to 10 x 10. Squares are QMC results for the ordinary single-band Hubbard 
models. The gray (dotted) curves show the Hartree-Fock results for our model (for the 
ordinary single-band Hubbard model). 
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Figure 2: Ground state energy per site at half filling as function of the oxygen distortion 
u. We take U/t = 4, g/t = 10 and k/t = 100. The curve is a polynomial fit to the 
extrapolated data to the limit iVs = oo. 
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Figure 3: The U dependence of (a) Jahn- Teller distortion u*, for the ground state at the 
optimized value of u — u*, (b) charge gap A c and (c) stabilization energy E JT . In (b), 
the charge gap is plotted for electron doping (A+) and hole doping (A~). Dashed lines 
show comparisons with experimental indications. 
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Figure 4: Inverse plot of peak values of the equal-time structure factor for the orbital 
correlation at the wavenumber (ir, ir) as function of doping concentration at U/t = 4. 
The gray line is the linear fit to the data at < 5 < 0.15. 
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Figure 5: Doping concentration dependence of chemical potential /j, at U/to — 4. The 
chemical potential is scaled by 5 2 . 
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Figure 6: Optical conductivity of the Mn model at U/to — 16 in the Lanczos diagonaliza- 
tion result of 14 electrons on the 4 by 4 lattice. The ground states are chosen by taking 
optimized phase shift of the boundary condition in x and y directions. The dashed curve 
D represents the Drude part. For the convenience of illustration, the delta functions ob- 
tained in the Lanczos diagonalizations are replaced with the Lorenzians with the width 
O.lto- 
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